#include <buola/mat.h>
#include <buola/mat/prettyprint.h>
#include <buola/mat/svd.h>
#include <buola/mat/traits_typename.h>

using namespace buola;

int main(int /*argc*/,char** argv)
{
    mat::CMat_d a=mat::random(2000,2000,std::normal_distribution<double>());
    mat::CMat_d b=a;
    //a(1500,1500)=0;

    start_timer();
    bool lEqual=true;
    for(int r=0;r<a.Rows()&&lEqual;r++)
    {
        for(int c=0;c<a.Cols();c++)
        {
            if(a(r,c)!=b(r,c))
            {
                lEqual=false;
                break;
            }
        }
    }
    end_timer();
    msg_info() << lEqual << "\n";

    start_timer();
    if(a==b)
        lEqual=true;
    else
        lEqual=false;
    end_timer();
    msg_info() << lEqual << "\n";
}
